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Abstract 



We present Herschel/PACS spectroscopic maps of the dwarf galaxy NC4214 observed in 6 far infrared fine-structure lines: [C ii] 
158;um, [O III] SSyum, [O i] 63yum, [O i] 146yum, [N ii] 122yum, and [N ii] 205yum. The maps are sampled to the full telescope spatial 
resolution and reveal unprecedented detail on ~ 150 pc size scales. We detect [C n] emission over the whole mapped area, [O iii] 
being the most luminous FIR line. The ratio of [O m]/[C ii] peaks at about 2 toward the sites of massive star formation, higher than 
ratios seen in dusty starburst galaxies. The [C ii]/CO ratios are 20 000 to 70 000 toward the 2 massive clusters, which are at least an 
order of magnitude larger than spiral or dusty starbursts, and cannot be reconciled with single-slab PDR models. Toward the 2 massive 
star-forming regions, we find that L[cii] is 0.5 to 0.8% of the Ltir. All of the lines together contribute up to 2% of Ltir. These extreme 
findings are a consequence of the lower metallicity and young, massive-star formation commonly found in dwarf galaxies. These 
conditions promote large-scale photodissociation into the molecular reservoir, which is evident in the FIR line ratios. This illustrates 
the necessity to move to multiphase models applicable to star-forming clusters or galaxies as a whole. 
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1. Introduction 

As building blocks of larger galaxies, dwarf galaxies in the early 
universe presumably played an important role in the evolution of 
galaxies we see today. Understanding how stars form and evolve 
under low metallicity conditions and the process of subsequent 
enrichment of the interstellar medium (ISM) will help us con- 
strain possible scenarios of galaxy formation and evolution. The 
Local Group dwarf galaxies provide laboratories to study in de- 
tail the effects of the lower metal abundance on the process of 
star formation and the feedback on the ISM. 

To examine the interplay between massive star formation 
and low metalUcity ISM, we have mapped the far infrared 
(FIR) fine structure lines in the Local Group Magellanic-type 
irregular galaxy, NGC4214, as part of the Herschel key pro- 
posal, SHINING (P.I. E. Sturm). NGC4214 harbors 2 main 
star-forming complexes containing hundreds of O stars as well 
as a super star cluster (SSC). Our aim here is to charac- 
terise the physical conditions and the structure of the different 
phases of the ISM in NGC4214. The proximity of NGC4214 
(2.9Mpc, iMaiz-A pellaniz et al. 200 ^ and its low -metallicity 
(Iog(0/H) + 12 ^8.2. iKobuInickv & SkiIIman|[T996 ) make it a 
choice target. The PACS spectrometer onboard the Herschel 
Space Observatory allows us to zoom into the low-metallicity 
photodissociation regions (PDRs) and their surroundings, trac- 
ing physical properties from 125 to 165 pc size scales. 

The tools for this study are the FIR fine-structure lines aris- 
ing from some of the most abundant species C, O, and N, which 



provide the most important cooling channels in the neutral and 
diffuse ionized ISM of galaxies. These lines trace H ii re- 
gions and PDRs, which are neutral atomic and molecular regions 
where penetrating FUV photons (6 eV < hv < 13.6 eV) domi- 
nate the energy balance and the chemical composition of the gas. 
Most of the ISM in galaxies is cont ained in PDRs. The FIR lines 
in combination with PDR models ("Tielens & Hollenbach| 
Wolfireetal. 1990; Ferland et al. 1998; Kaufman et al 



1985; 



19991: 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Abe l et al.1 120051: iRoIIig et al . 2006) can unveil the structure of 
these regions, which depends on both the strength of the FUV 
radiation field density (Go) illuminating the cloud and the hy- 
drogen density (nH). 



2. Observations and data reduction 

We mapped the 6 FIR fine-structure lines, [C ii] 158, [O iii] 88, 
[O i] 63, [O i] 146, [N ii 1 122 and [N ii] 205yum, using th e PACS 
spectrometer (Poglitsch et alJl2010l) on Herschel (Pilbra tt et al.l 
1201 Oh . The PACS' array consists of 5 x 5 spatial pixels of 9.4" 
size offering a total field-of-view of 47" on the sky. Each spa- 
tial pixel covers 16 spectral elements with a spectral resolution 
ranging from 100 to 300km/s. The lines are not spectrally re- 
solved for NGC4214. The 6 lines were observed in raster map- 
ping mode over a total of 15h during the Science Demonstration 
Phase. The [O iii] 88;um and [O i] 6'ipva ("blue") were mapped 
with 5x5 rasters and the other lines ("red") with 3x3 rasters, all 
centered on 12hl5m39.08s,+36dl9m35.9s (J2000). Each raster 
pointing is separated by approximately half of the array (22", 
24") in the "red" range and approximately one third of the array 
(14.5", 16") in the "blue" range, giving homogeneous coverage 
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Table 1. FIR line ratios 



NGC42I4 



total 



I 



II 



III 



Figure 1. Individual PACS spectral lines added up over the cen- 
tral raster The lines are not spectrally resolved. 

of the galaxy and a total field of view of 1.6'xl.6'. We selected 
the chop-nod mode with a chop throw of 6' off the source. 

The data were reduced with the PACS spectrometer pipeline 
of the Herschel Interactive Processing Environment (HIPE) 
v3. 0.455. We show the 6 individual lines in Fig. [1] We fitted 
a spline plus Gaussian to the line and baseline and measured the 
Icr noise of the residual. In this way, we estimate the overall 
noise in the maps to be 0.17, 0.42, 0.27, 0.07, 0.06, and 0.30 
Jy/pix for the [C ii] 158yum, [O iii] 88//m, [O i] 63yum, [O i] 
146yum, [N ii] 122/vm, and [N ii] 205jum lines, respectively. We 
applied a flat field correction and divided by calibration correc- 
tion factors of 1.1 ("red") and 1.3 ("blue") to obtain the final 
reduced maps. 

The rasters are combined by a drizzle scheme to produce fi- 
nal maps with 4" grid pixel size for the [O iii] and [O i] 63/vm 
lines, and 6" for the other lines. To compare the lines and create 
ratio maps, we scaled all maps down to the coarser resolution of 
6". We analyze separately the whole map as well as three spe- 
cific regions corresponding to NGC4214-I (center), NGC4214- 
II (southeast), and NGC4214-III (northwest). These are repre- 
sented by the red circles on Fig.|2^. Regions I and II are centered 
on the [C ii] peaks, with radii encompassing most of the emisson 
(at the level of ~50% of the peak), while region III was chosen 
by taking into account the peak of the CO emission as well as the 
brightest counterpart of the [C ii] emission (S/N > 15), which 
is truncated by the map edge. 

3. Results 

The HerschelfPACS maps clearly resolve the two star-forming 
regions NGC4214-I and NGC42 14-11, where most of the emis- 
sion originates. The [N ii] 122/im emission is rather faint and the 
[N ii] 205//m is not detected. The most prominent lines are the 
[C ii] 158//m, [O III] 88//m, [O i] 63jum, and, to a lesser extent. 



RA[J2000] 1215 39.1 1215 39.1 1215 40.5 12 15 36.7 

Dec[J2000] +3619 35.9 +3619 35.9 +3619 05.9 +362001.0 
Aperture (") - 18 12 9 



Quantity 



Line ratios 



/[0IlIj//[ClI] 


0.996 


1.40 


1.71 


0.264 




(0.009) 


(0.007) 


(0.014) 


(0.053) 


/[0I]J146//[0iy63 


0.058 


0.060 


0.072 


- 




(0.015) 


(0.011) 


(0.023) 


(0.095) 


/[0I]i63//[CIll 


0.445 


0.474 


0.560 


0.438 




(0.007) 


(0.004) 


(0.008) 


(0.035) 


/[0I]J146//[CII] 


0.026 


0.028 


0.040 


- 




(0.007) 


(0.005) 


(0.013) 


(0.042) 


Ltir [10' Lq] 


37.7 


8.60 


4.52 


0.522 


10-3(L[oiy(,3 


2.77 


8.11 


6.50 


11.1 


+ -'^[CII]) / -^TIR 










10-3L[cii]/Ltir 


5.30 


7.93 


5.71 


7.70 


10 ^[cii]/^co 


33.8 


75.0 


22.9 


4.73 



the [O i] 146/ym and [N ii] 122yum lines. The emission from the 
[C n] line is the most extended of all of the lines (Fig. |2^) and 
covers the entire map. The total [C ii] luminosity in the PACS 
map is 2.0 X 10^ Lq with about 50% of the emission originating 
in the two bright star-forming sites and 50% in the surrounding 
medium. The [C ii] 158/im line is known to be one of the most 
important tracers of physical conditions in PDRs because it is so 
luminous and has a low critical density for collisional excitation. 

Comparing the [C ii] and [O i] 63/im line emission, we 
find that the latter is a factor of two weaker integrated over the 
mapped region. These two lines show a very similar distribution 
(Fig.|2j5), both peaking at the 2 star-forming sites and spanning 
a small range of ratios. These ratios vary by less than a factor 
of two across the entire map. The [C ii] line, in principle, may 
originate in the PDR and the diff'use media alike. The observed 
resemblance between the [O i] 63;um and [C ii] maps (Fig. |2j5) 
suggests that the [C ii] is dominated by the PDR component (see 
also Sect. [TB. 

The [O i] 146;um line is fainter than the [O i] 63yum line, and 
is barely detected in the outer regions of the map. Toward the 
2 peaks, and in total, the ratios of the two [O i] lines are about 
0.06 (Table [TJ, and do not vary much across the map. These ra- 
tios are insensitive to density since their critical densities are not 
very different. For moderate densities (n < 10'' cm"^), this ratio 
is more indicative of the gas temperature. 

The [O III] 88/im line is the brightest of all lines observed in 
NGC4214 (Fig. |2j;). A comparison of the [O m] line map with 
the HS T optical images of NGC4214 (lUbeda et al.ll2007h . illus- 
trates that it coincides with the star formation sites and exhibits 
a peak toward the position of the SSC in NGC4214-I, thus trac- 
ing the ionization source. High ratios of [O m] to [C ii] observed 
here (Table[T]i are not commonly seen in dusty starbursts or spiral 
galaxies ( Negishi et al...2001: .Malhotra et al...2001.) . an example 
being M82 with a ratio of ~ 0.7. However, high ratios have been 
observed in the ISO observations of dwarf galaxies and may 
be explained by ho t stars surrounded by optically thick PDRs 
([Hunter etal.ll200lh . The [O in] / [C ii] ratio map of NGC4214 
(Fig. |2]l) exhibits peaks in intensity toward the 2 star-forming 
sites but drops off rapidly away from the peaks. We note that the 
peaks of the [O iii] / [C ii] ratio map appear to be offset from 
those of the [C ii] by ~ 6", which is larger than the expected rel- 
ative pointing uncertainty (on the order of several arcsecs). The 
shift in the [O iii] / [C ii] peak, corresponding to a shift of ~ 80 
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Figure 2. PACS spectral maps of NGC4214. Fig 2a: PACS map of the [C ii] 158yum line with CO(l-O) contours (IWalter et al.1 
1200 ih . The total field of view is 1.6'xl.6' with a grid resolution of 6" . The red circles show the positions and sizes of the 3 regions 
we study in more detail: NGC4214-I (center- 36" diameter), NGC4214-II (southeast-24"), and NGC4214-III (northwest- 18"). The 
dotted circle on the lower right corner represents the PACS beam size of 9.4"x9.1" at 62yum, 9.6"x 8.8" at 90jum, and 11.8"x 11.0" 
at 154/im. Fig 2b: PACS map of the [O i] 63jum line with [C ii] 158jum contours. Fig 2c: PACS map of the [O iii] 88jum line with 
[C ii] 158;um contours. Fig 2d: Ratio map of [O iii] 88yL(m/[C n] with [C ii] contours. Units for the color bars are 10"'^ W vaT^ grid 
pix"'. The [C ii] contours begin at 0.17 x 10"'^ W m"^ (above 15cr) and are linearly spaced up to the peak of 1.66 x 10""' W vnT^. 



pc, if real, implies that the ratio peak coincides with the SSC 
peak, the probable source of the ionization. Since O^^ has an 
ionization potential of 35.1 eV, [O in] traces the highly ionized 
medium around early O stars and is thus more confined than the 
[C n] emission. This is consistent with the AKARI findings that 
the [O III] 88/im line traces excitation sources even more accu- 
rately than the radio continuum (Okada et al. 2009) . 



We show the CO contours from IWalter et aP (1200 ll) in 
Fig. |2^. There is a striking distinction between the three dif- 
ferent regions in terms of [C ii]/CO (Table [D- NGC4214-III 
probably hosts little star formation, which explains its paucity in 
excited FIR line emission compared to its molecular reservoir. 
However, even the two main regions of active star formation ex- 
hibit very different [C n]/CO ratios. One explanation of these 
high [C ii]/CO ratios is the presence of small clumps of CO em- 
bedded in large PDR envelopes. Applying this interpretation to 
these observations would imply that the star formation in the 
central region of NGC42 14-1 has (already) affected (photodisso- 
ciated) most of the surrounding molecular reservoir, while many 
of the molecular clouds in the south eastern region are relatively 
unscathed. It is interesting to remark that NGC42 14-11 has a 
younger stellar population than NGC4214-I dUbeda et al.ll2007h . 



Using Spitzer data, we estimate the total infrared lumi- 
nosity (Ltir) over the area mapped with PACS, following 
iDale & Heloul (l2002h (see Table Q}. The [C n] Une alone ac- 
counts for ~ 0.5% of the Ljiio which is at the high end of 
the range of normal galaxi es and starburst galaxies (IStacev et al.l 
I1991I: iNegishi et al.ll200ll) and more typical of the elevated ra- 
tios seen in dwarf galaxies ( Madden 2000 ; Hunter et al. 200 ij ). 
Toward the two bright star-f orming regions, Lfcni / LxiR is sim- 
ilar to 30 Dor in the LMC (iPoglitsch et al.lll995i- iRubin et al.1 
|2009|) . Likewise, L[oiii]/iTiR either toward these 2 prominent 
star- forming regions or over the entire galaxy is an order of mag- 
nitude higher than what'Ne gishi et al.) (i2001 ) find for normal and 
starburst galaxies. The FIR lines of the 2 star-forming regions 
contribute all together up to about 2% of their TIR flux. 



3.1. Origin of ttie [C n] line 

Before we can consider PDR model solutions to the observa- 
tions, we consider the origin of the C^ line, which is most often 
used as a diagnostic of PDR conditions because it is a primary 



' To compare rat ios using Lfir in the literature, for dwarf galaxies 
jHunter et alj|200lh . Ltir is approximately a factor of 2 greater than 
Lfir, as defined by IRAS 60 and lOOyum fluxes dHelou et alil988h . 
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coolant in PDRs. However, since the ionization potential of car- 
bon (11.3 eV) is less than that of hydrogen, C^ can exist in both 
the ionized and neutral gas. In ga lactic nuclei, most of the [C n] 
emission arises from the PDRs ( Stacev et al.ll99ll : lNegishi et alJ 
12001 ). With an ionization potential of 14.53 eV, the N^ emission 
originates only from the H ii regions. Therefore we use the [N ii] 
line to coiTect for a possible contribution of the ionized gas to 
the C+ Une. 

The [N ii] 205;um line was not confidently detected and could 
not be analyzed but the [N ii] I22jj.m line is detected in the two 
star-forming regions. We see evidence of [N ii] llljim emis- 
sion over the central 3x3 pixels in each region and average them 
to obtain a higher signal-to-noise ratio and measure the inten- 
sity of the line in these two regions. The observed ratios of the 
[N II] 122|/m/[C II] 158yum are 0.04 and 0. 03 in NGC42 14-1 and 
NGC42 14-11, respectively. Following Mal hotra et"an (1200 ih . we 
estimate a maximum contribution of the ionized medium to the 
[C ii] emission to be a factor of 1 1 times the flux of the [N ii] 
122yum. Therefore, at least 60% of the [C ii] emission originates 
in PDRs. To bound the physical conditions in NGC4214, we 
model the PDR (see Sect. 13.2b assuming that all [C ii] emission 
is from PDRs and that only 60% originates in PDRs. 

3.2. Results from PDR models 

As a first diagnostic of PDR properties, we use the model of 
iKaufman et al.l d 19991) by means of the web interface of the PDR 
Toolbojqj. This model, which assumes solar metallicity, calcu- 
lates the intensity of the main PDR lines as a function of the 
physical conditions, which are the incident radiation field (Go) 
and the density («h)- We can compare our observed line inten- 
sities to the model predictions to narrow down the range of the 
physical conditions. 

We use [O i] 63//m/[C ii], [O i] 146;um/[C ii], 
[O i] 146/im/[0 i] 63;um, and the ([C ii]h-[0 i] 63//m)/lTiR 
and [C ii]/CO (Table [B. The ratio ([C ii]-h[0 i] 63jum)/l7-/R is 
interesting because it relates the "total" cooling of the gas to the 
total cooling by the dust, so measures the efficiency of the dust 
to heat the gas (photoelectric efficiency). The observed [C ii]/CO 
ratio disagrees significantly with those predicted by the best-fit 
solution for all the other line ratios. This has been noted before 
for other low metallicity PDRs (see Sect. 13.31 1 and may be a 
geometry or low metallicity effect not taken into account in 
the model. Using all the above line ratios except [C ii]/CO, the 
maximum values found for NGC4214-I are Gcu ~ 800 and 
nn ~ 2000. For NGC4214-II, we find that Go ~ 1000 and 
wh ~ 3000. The solution for the third region NGC4214-III infers 
a lower incident radiation field (Go < 100). The quoted values 
above take into account the contribution to the [C ii] intensity of 
the ionized medium between and 40%. In all cases, these PDR 
parameters are indicative of a relatively moderate gas density 
and radiation field. 



3.3. [C ii]-to-CO ratio 

We measure [C ii]/CO intensity ratios of 70000 and 20000 for 
NGC4214-I and II, respectively (Table [TJ. These values are at 
least an order of ma gnitude higher thaii those inferred for sp iral 
or starburst galaxies lstacev et al.lll99 1':'Negishi et al."2001') but 
in clo ser agr eement with those of dwar f galaxies (Poglitsch et al. 
Il995t iMadden et alJll997t iHunter et al.ll2001h . NGC4214-I has 



http://dustem.astro.umd.edu/pdrt/index.html 
in units of Habing field: 1.6 x 10"' erg cm"- s" 



a [C i i]/CO similar to 30 Dor (iPoglitsch et al.lll995HRu"bin et al.l 
12009 =). This high ratio could be explained by the paucity of dust; 
less dust shielding allows the photodissociating FUV photons 
to penetrate deeper into the cloud, reducing the size of the CO 
cores while leaving a large [C ii] -emitting envelope. This en- 
velope may well contain self-shielded molecular hydrogen. If 
this were true, the [C ii] would be a more accurate tracer of the 
molecular reservoirs in low metallicity environments than CO. 

4. Conclusion 

These observations of NGC4214 with Herschel/PACS have en- 
abled us to spatially resolve the FIR fine-structure lines in the 
low metallicity dwarf galaxy NGC4214. A very striking re- 
sult from the mapping is the brightness of the [C ii] 158jum , 
[O i] 63/zm, and the [O iii] 88//m lines. The [C ii] alone cor- 
responds to 0.5% of the total infrared luminosity. The com- 
bined luminosity of the FIR lines reaches 2% of the Ltir. The 
[O III] 88/im line is more confined to star-forming regions and 
the [O iii]/[C ii] reaches almost a factor of 2 close to the SSC in 
NGC4214-I. 

The [C ii]/CO ratios are also exceptionally high (20000 - 
70000). These extreme ratios are indicative of an ISM that is 
strongly affected by photodissociation and reflect the combined 
effects of low metallicity ISM and intense star formation. It 
should be borne in mind that, when comparing with models, 
the models predict line ratios for a single, homogeneous PDR. 
In reality, we observe a multiphase medium that probably har- 
bors many extremely different PDR conditions. Therefore, this 
preliminary comparison indicates the need for more complete 
(complex) models to be able to derive accurate physical quanti- 
ties. 
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